Dilepton production by dynamical quasiparticles in the strongly interacting 

quark gluon plasma 
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The dilepton production by the constituents of the strongly interacting quark-gluon-plasma 
(sQGP) is addressed. In order to make quantitative predictions at realistically low plasma tem- 
peratures (0(T C )), experimentally relevant low dilepton mass (0(1 GeV)) and strong coupling 
(as ~ 0.5-^1), we take into account not only the higher order pQCD reaction mechanisms, but 
also the non-perturbative spectral functions and self-energies of the quarks, anti-quarks and gluons 
thus going beyond the leading twist. For this purpose, our calculations utilize parametrizations of 
the non-perturbative propagators for quarks and gluons provided by the dynamical quasi-particle 
model (DQPM) matched to reproduce lattice data. The DQPM describes QCD properties in terms 
of single-particle Green's functions (in the sense of a two-particle irreducible approach) and leads 
to the notion of the constituents of the sQGP being effective quasiparticles, which are massive 
and have broad spectral functions (due to large interaction rates). In the present work, we derive 
the off-shell cross sections of dilepton production in the reactions q + q — > l + l~ (Drell-Yan mech- 
anism), q + q — > g + l + l~ (quark annihilation with the gluon Bremsstrahlung in the final state), 
q{q) + 9 q{q) + ' + ' _ (gluon Compton scattering), g — > q + q + l + l~ and q(q) — > q(q) + g + l + l~ 
(virtual gluon decay, virtual quark decay) in the sQGP by dressing the quark and gluon lines in the 
perturbative diagrams with the DQPM propagators for quarks and gluons. 

PACS numbers: 25.75.Nq, 25.75.Cj, 24.85.+p, 12.38.Lg, 13.60.Hb 



I. INTRODUCTION 



Since many years the transition between the hadronic 
phase and the quark-gluon plasma (QGP) as well as 
the nonperturbative properties of QGP motivate a large 
community and justify large-scale experiments, in which 
heavy nuclei are collided at relativistic energies in or- 
der to achieve the high energy-densities necessary for the 
transition to the deconfined state of matter. Electromag- 
netic probes (i.e. dileptons and photons) are powerful 
tools to explore the early (hot and dense) stage of the 
heavy-ion collision, since, unaffected by the final state 
interaction, they carry to the detector information about 
the conditions and properties of the environment at the 
time of their production - encoded in their mass and mo- 
mentum distributions, - thus providing a glimpse deep 
into the bulk of the strongly interacting matter (3, @| . In 
1978, E. Shuryak proposed to use dilepton as probes of 
QGP, after the suggestion was made that the dileptons 
and photon yields reflect the properties of the medium 
created in hadron-hadron collisions (see the pioneering 
works 0-0|). 

Real and virtual photons, i.e. dileptons, are emitted 
over the entire space-time evolution of the heavy-ion col- 
lision, from the initial nucleon-nucleon collisions through 
the hot and dense phase and to the hadron decays af- 
ter freeze-out. This is both a challenge and advantage of 
electromagnetic probes. Fortunately, lepton pairs possess 
an additional degree of freedom (the invariant mass Q 2 ), 
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which allows to separate different "physics" by observing 
the dilepton radiation in different mass ranges. The low 
mass (Q < 1 GeV) spectrum of dileptons - generated 
in heavy ion collisions - is dominated by the vector me- 
son decays, the production of lepton pairs of high mass 
(Q > 3 GeV) is governed by the perturbative quantum 
chromodynamics (pQCD), while the dilepton yield in the 
intermediate mass range (1 < Q < 3 GeV) is sensitive to 
the possible formation of a QGP. 

Dilepton measurements have possibly provided a sign 
of the deconfined matter at SPS energies. The NA60 
Collaboration [8T4l0| has recently found that the effective 
temperature of the dileptons in the intermediate mass 
range is lower than the T e g- of dileptons at lower masses, 
which arc of hadronic origin. This can be explained, if 
one assumes that the spectrum at the invariant masses 
above 1 GeV is dominated by the partonic channels in the 
QGP [ITI4T3I ] . In this case, the softening of the transverse 
mass spectrum with growing invariant mass implies that 
the partonic channels occur dominantly before the col- 
lective radial flow has developed. The assumption that 
the dilepton spectrum at masses above 1 GeV could be 
dominated by QGP radiation was supported by the stud- 
ies within the Hadron-String-Dynamics (HSD) transport 
approach (l4^ . which have shown [l5| that the measured 
dilepton yield at low masses (Q < 1 GeV) can be ex- 
plained by the dilepton production in the hadron inter- 
action and decay, while there is a discrepancy between 
the HSD results and the data in the mass range above 
1 GeV. This access seen at Q > 1 GeV is not accounted 
for by hadronic sources in HSD - in-medium or free - 
and might be seen as a signal of partonic matter, mani- 
fest already at 158 AGeV incident energy. 



(a) (6) (c) (rf) 

FIG. 1: Perturbative QCD diagrams contributing to the dilepton production up to the order O(ag): (a) Drell-Yan mechanism, 
(b) gluon Compton scattering (GCS), (c) vertex correction, (d) quark annihilation with gluon Bremsstrahlung. Virtual photons 
(wavy lines) split into lepton pairs, spiral lines denote gluons, arrows denote quarks. In each diagram, the time runs from left 
to right. 



Recently, the PHENIX Collaboration has presented 
first dilepton data from pp and Au + Au collisions at 
Relativistic-Heavy-Ion-Collider (RHIC) energies of ^/s = 
200 GeV [l6l - fl9l |. The data show an even larger en- 
hancement over hadronic sources in Au + Au reactions in 
the invariant mass regime from 0.15 to 0.6 GeV, which 
could not be explained in the scope of the HSD approach 
neither by meson decays - in-medium or free - nor by 
hadronic Bremsstrahlung [lj|. It is of interest, whether 
the excess at RHIC is due to the dominance of sources 
in the QGP [Hj. 

Early predictions of the dilepton emission from QGP 
relied on perturbative formulae for the cross sections of 
the virtual pho ton production in q + q and q + g colli- 
sions {3. l20l. l2ll| . Indeed, first concepts of the QGP were 
guided by the idea of a system of partons which interact 
weakly, with pQCD cross sections. However, most theo- 
retical estimates of the temperatures, which are reason- 
ably expected to be currently achieved in heavy ion colli- 
sions are not extremely large compared to the QCD scale 
Aqcd 22]. Therefore, the QCD coupling as is not small. 
In agreement with this early expectation, experimental 
observations at RHIC indicated that the new medium 
created in ultra-relativistic Au+Au collisions was inter- 
acting strongly - stronger than hadronic matter. More- 
over, in line with theoretical studies in Refs. [23T - [25| the 
medium showed phenomena of an almost perfect liquid 
of partons [2q - l29| as extracted from the strong radial 
expansion and elliptic flow of hadrons [26j-|29(. Studies 
performed in the framework of the lattice regularized 
QCD [30J have also shown that the high temperature 
plasma phase is a medium of interacting partons which 
are strongly screened and influenced by non-perturbative 
effects even at temperatures as high as 10 T c . 

Consequently, the concept of perturbatively interact- 
ing quarks and gluons as constituents of the QGP had 
to be given up. Due to large running coupling, the next- 
to-leading order (NLO) gluon-quark interactions are ex- 
pected to contribute considerably in addition to the lead- 
ing order mechanism of quark-quark annihilation (qq — ¥ 



l + l ) to the QGP radiation spectrum. The importance 
of the higher order corrections is long understood 
On the other hand, non-perturbative nature of the sQGP 
constituents manifests itself not only in their strong cou- 
pling, but also in the modified spectral densities and self 
energies, which should be taken into account in consistent 
calculations of dilepton production from the QGP. 

A cure can be found in reordering perturbation the- 
ory: by expanding correlation functions in terms of effec- 
tive propagators and vertices instead of bare ones [32l |. 
A powerful resummation technique was developed by 
Braaten, Pisarsky [33[ and Wong [34]]. The production 
of dileptons was calculated at leading order in the effec- 
tive perturbation expansion in (33|, using as the effec- 
tive propagators the bare ones plus one loop corrections 
evaluated in the high-temperature limit [35j-|38[. In this 
approach the singularity of the production cross section 
- that dominates the dilepton rate - is regularized by 
the thermal masses of quarks mth and gluons m g , which 
are in turn determined by the one-loop leading order re- 
sult in the thermal perturbation theory (HTL). The ap- 
proach has been extended to the dilepton radiation from 
non-equillibrium plasmas in [39l |40| . 

However, since virtual photon rates need to be evalu- 
ated at temperatures that are not very large compared 
to T c , it is important to take values for m t h, tn g not 
from the HTL approximation but from, for instance, a 
fit of the lattice QCD entropy by a gas of massive quarks 
equation of state 0, HU , as has been done in [43|, 13 ■ 
Alternatively, one might treat thermal masses in the cal- 
culation of the dilepton rates as phenomenological pa- 
rameters [4a |. 

Not long ago, a first attempt appeared to calculate 
directly on the lattice the production of dileptons in 
QGP [46j . The suppression at small Q 2 observed on the 
lattice has attracted a lot of interest, because it is not 
what one would expect from (resummed) perturbation 
theory: The finite thermal masses would indeed produce 
a drop of the Born (Drell-Yan) term Q+q — > 7* because of 
the threshold effect - as predicted 47J in relation to the 



cut-off in the momentum distribution of quarks and con- 
firmed in effective perturbation theory in the works [43l — 
Hil, - but there are higher order processes (qq — > j*g, 
19 ~~ * 7*9) that have no threshold and would fill the 
spectrum at small Q 2 . The effect of multiple scattering 
of the quark in the plasma (Landau- Pomeranchuk-Migdal 
effect [48 50]) on the rate of qq — > 7* was considered in 
the context of the semi-classical approximation in [5l| . 
where it was shown that the multiple scatterings lead to 
the disappearance of the threshold behavior in the Born 
rate. 

In the present work, we pursue a model, in which the 
(multiple) interaction of quarks and gluon in sQGP is 
encoded in their effective broad spectral functions. The 
non-zero width of quarks is related to their strong inter- 
action, which is manifest in the elastic scattering as well 
as in the virtual gluon emission. By dressing the partonic 
lines in the dilepton production processes (q + q — > l + l~ , 
q + q -+ g + l + l~ , q + g -+ q + l + l~) with effective spectral 
functions we study the effect of the partonic interactions 
in the plasma on their dilepton radiation, especially in 
the interesting region of low Q 2 . 

The effective propagators were obtained from lattice 
data in the scope of the Dynamical QuasiParticle Model 
(DQPM) [H. DQPM describes QCD properties in terms 
of single-particle Green's functions (in the sense of a two- 
particle irreducible (2PI) approach) and leads a quasi- 
particle equation of state, which reproduced the QCD 
equation of state extracted from Lattice QCD calcula- 
tions. According to the DQPM, the constituents of the 
sQGP are effective strongly interacting massive partonic 
quasi-particles with broad spectral functions, i.e. non- 
vanishing width. 

Strong interaction of partons - reflected in their broad 
widths - leads to higher-twist corrections to the standard 
pQCD cross sections [53j. The influence of the higher 
twists is hard to estimate and has been ignored in most 
calculations so far. Limiting the calculations to the lead- 
ing twist is a widely used approximation at high energies 
for the following reason: The higher twists by definition 
are vanishin g in the limit of infinite invariant mass of the 
lepton pair [5J]. Indeed, one observes that by calculat- 
ing next-to-leading order contributions and refitting the 
parton distributions accordingly [3ll . |55| , one can already 
significantly reduce the discrepancy with the data on the 
double-differential Drell-Yan cross section in elementary 
p + p collisions [56l - l58j compared to the LO predictions. 

However, the power-Q 2 suppressed (i.e. higher twist) 
contributions can be large at realistic energies and/or 
in collective systems with a relatively short interaction 
length like the sQGP created in heavy-ion collisions. It 
has been shown [53J] that the contribution of higher twists 
is essential for a proper description of the data on the 
triple- differential Drell-Yan cross section. For instance, 
quark and gluon off-shcllncss - arising due to the non- 
perturbative interaction between the partons - have a 
large effect on the transverse momentum distribution of 
produced lepton pairs [H, [H, H3| • 
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FIG. 2: Diagrams contributing to the dilepton production 
by virtual quasi-particles in addition to the ones presented 
in Fig. [JJ Upper part: the decay of a virtual quark; Lower 
part: the decay of a virtual gluon. Virtual photons (wavy 
lines) split into lepton pairs, spiral lines denote gluons, arrows 
denote quarks. 



In the current paper, we study the dilepton production 
from dynamical quasi-particle interactions. For this pur- 
pose we derive the off-shell cross sections of q + q -+ l + l~, 
q + q^ g + l + l~, q + g -> q + (q + g -> q + l + l~), 
Q -> q + g + l + l~ (q -> q + g + l + l~) and g -> q + q + l + l~ 
by calculating them for the arbitrary virtualities of exter- 
nal quarks and gluons, while dressing the internal lines 
with effective self energies, and take into account the 
non-zero widths of partons by convoluting the obtained 
cross sections with the effective spectral functions from 
the DQPM. Using the DQPM parametrizations for the 
quark(gluon) self-energies, spectral functions and inter- 
action strength, we calculate the dilepton production by 
partonic interactions in the sQGP while accounting for 
their non-perturbative dynamics, including the higher 
twist effects. 

The perturbative diagrams for the dilepton production 
at order up to 0(as) are illustrated in Fig. [TJ Let us 
briefly summarize the differences of our effective theory 
approach from the standard pQCD: 

• (a) we take into account full off-shell kinematics, in 
particular the transverse motion and virtuality of 
the partons, 

• (b) quark and gluon lines are dressed with non- 
perturbative spectral functions and self-energies: 
the cross sections are derived for arbitrary masses 
of all external parton lines and integrated over 
these virtualities weighted with spectral functions 
(see e.g. Refs [53|, E[ f° r an introduction to the 
method); the internal lines are dressed with self 
energies, 

• (c) vertices are modified compared to pQCD by re- 
placing the perturbative coupling (that runs with 
the momentum transfer) with the full running cou- 
pling as that depends on the temperature T ac- 
cording to lattice data parametrization by [52j . 
while the temperature is related to the local energy 
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density e by the 1QCD equation of state. Note that 
close to T c the full coupling increases with the de- 
creasing temperature much faster than the pQCD 
prediction. 

• (d) Due to the broad width for quarks and glu- 
ons in the sQGP [43[ - which is the consequence 
of their high interaction rate,- we obtain non- 
vanishing contributions also from the processes of 
the decays of virtual quarks (q —> q + g + l + l~) 
and gluons (g — > q + q + l + l~), which are forbidden 
kinematically in pQCD (see Fig. [2]). 

Thus, we utilize the DQPM parametrization for the ef- 
fective quark and gluon propagator while simultaneously 
taking into account the NLO reaction mechanisms. This 
way, we are going in our calculations beyond the leading 
order in as and depart from the leading twist approxi- 
mation. 

The paper is organized as follows. The analysis of the 
off-shell kinematics and the calculations of the off-shell 
cross sections in our model are given in section IIIII for 
each of the partonic processes separately. On the other 
hand, in the limit of high hard scale Q 2 , the off-shell 
cross sections should approach the perturbative ones 
l62j . Therefore, let us first recapitulate the corresponding 
pQCD results in section HU The off-shell cross sections 
will be compared to the perturbative ones throughout 
the section iDll 

In Section IIVI we give a simple example of an applica- 
tion of the off-shell cross sections derived above to the 
calculation of the relative contributions of different pro- 
cesses to the sQGP radiation. We point out, however, 
that the considerations in the current paper will prob- 
ably be not detailed enough to be applied to the de- 
scription of the hot and dense early phase of relativistic 
heavy-ion collisions, where sQGP is formed. A quan- 
titative comparison to the experimental data and reli- 
able conclusions on the relative contribution of various 
sources to the experimentally observed thermal dilepton 
spectrum (§4l0l [l6l - [l9j requires taking into account the 
non-equilibrium dynamics of the heavy-ion collision in its 
full complexity by use of microscopic transport models, 
which is beyond the scope of this study. The main pur- 
pose of the current paper is to built an effective approach 
for the derivation of the off-shell cross sections for the in- 
teraction of dynamical quasi-particles as constituents of 
the sQGP. The qualitative analysis of the relative im- 
portance of different processes in section IIVI should be 
understood as an illustration to the above results rather 
than a quantitative prediction for the dilepton yield from 
heavy-ion collisions. 

Section [V] is devoted to analyzing the effect of finite 
quark and gluon widths on the dilepton rate. In Sec- 
tion [VI] we summarize the main results and their possible 
applications. 



II. DILEPTONS FROM PERTURBATIVE 
PARTONS 

In the present Section, we consider the following par- 
tonic mechanisms for dilepton production in pQCD: 

1. the Drell-Yan mechanism of quark annihilation (q+ 

2. quark + anti-quark annihilation with gluon 
Brcmsstrahlung in the final state (q + q — > g + 7*), 

3. Gluon Compton scattering (q + g — > 7* + q and 
<7 + 9 -> 1* + <?)■ 



A. Drell-Yan mechanism 

The leading order pQCD mechanism for the dilepton 
production in the partonic phase is the same as for the 
well known Drell-Yan (DY) process (63[: quark and an- 
tiquark annihilate into a lepton pair {qq — > l + l~), as 
presented by the diagram (a) in Fig[T] The leading order 
leading twist (IT) pQCD result for the cross section of 
DY dilepton production is 

d 3 a{qg ^ l+l-) \ BY = 4 ™ 2e g x lX2 
dQHx F dq* T ) on _ sh ~ 9Q 4 x 1+ x 2 [ XlX2 > 

xS(q 2 T )S {Q 2 - x lX2 S NN ) S (x F - * 2 ~ Xl \) 

\ -L X\X 2 / 

where a is the electromagnetic fine structure constant, 
e q the fractional quark charge, the subscript "on-shell" 
stands for 'leading twist', the lepton pair has invariant 
mass Q 2 and transverse momentum qx- 

In collinear pQCD, the off-shellness, mass and trans- 
verse momentum of the annihilating quark and antiquark 
are neglected, and, therefore, the incoming parton mo- 
menta are related to the momenta of colliding nuclei as 
Pqtq) = xiPaIA. In this case, the parton momentum 
fractions x\ and xi are related to the virtuality and xf 
of the produced photon as (cf. delta- functions in (fT])): 

Q 2 = s = xix 2 S N n; (2) 

XF = (X2 - Xl)/(1 - Xl X 2 ). (3) 

Note that the denominator of the xf definition in ([3]) is 
omitted in some works, where an approximate definition 
x F ~ 2p z / v / S NN is used instead of x F = q z /{qz)max, 
s denotes the invariant energy for the partonic process, 
Snn ~ for the hadronic one. The kinematical limits for 
this process are 

Snn>Q 2 , \x f \ < l,s = Q 2 . (4) 
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B. Gluon Bremsstrahlung 



The cross section of the gluon Bremsstrahlung process 



qq 



g + h+li- is [6J, |6£ 



d 2 a(qq -> gl + l~ 
dQ 2 d cos 6 



&a e q as s 



Q 2 



-shell 



27Q 2 



1 + cos 2 9 + 4 



s 2 sin 2 6 
Q 2 s 



(s-Q 



2\2 



,(5) 



where s is the total energy squared of the colliding par- 
tons, and is the scattering angle of the outgoing lepton 
pair with respect to the incoming quark momentum in 
the quark center-of-mass system (CMS). Note that the 
cross section ([5]) can be written in terms of the Mandel- 
stam variables s, t and u as [20j 



d 2 a(qq -> gl + l~ 
dQ 2 dt 



shell 



8a 2 e 2 q a s [t - Q 2 ) 2 + (u - Q 2 ) 2 
27Q 2 s 2 tu 

xS(s + t + u-Q 2 ) (6) 

t u 2sQ 2 ' 
it t 

' (7) 



Sa 2 e 2 a s 



27Q 2 .s 2 V" ' * ' tu 
x 5(s + t + u-Q 2 



which coincides with the QED cross section for the vir- 
tual Compton scattering up to the color factor and the 
crossing transformation 66]. Here we denote the mo- 
menta of the incoming quark and antiquark as p\ and 
P2, the momenta of the outgoing gluon as virtual photon 



as k and q, s = (P1+P2) , t = (pi-qY 



(p2-q) 2 - The 



^-function S(s + u + t — Q 2 ) reflects the on-shell condition 
for the partons: 



Pi 



pi 



0. 



(8) 



The collinear divergence of the gluon Bremsstrahlung 
cross section for t —> and u — > (i.e. cos0 —> ±1) 
is obvious; a cut-off A 2 on t can be used in order to 
regularize it: 



t < 
t > 



-A 2 



-s + Q 2 +A 2 . 
Since in the CMS of the colliding partons we have 

t = Q 2 - + Vs\q\ cosQ, 
the corresponding cut-off with respect to cos is 

(V^q° -Q 2 - a 2 ) 



I cos 61 < 



(9) 
(10) 



(11) 



(12) 



Another divergence in the perturbative expression ([6]) is 
the infrared (IR) divergence for the energy of the gluon 
fc° — > due to the vanishing quark and gluon masses. 
Indeed, if all the partonic masses (virtualities) are ne- 
glected, we have in the CMS: 



t = k Q {-y/s~ + y/s~cosQ 2 ) -> at k° ->■ 0. 



(13) 



This divergence can be remedied by introducing a small 
finite gluon mass fi cut (cf. the plasmon mass in _22|). 
Indeed, the gluon thermal mass fi plays the role of a 
natural cut-off in the sQGP (cf. section UlI Bp . 



C. Gluon Compton scattering 

In QED, the Compton process refers to elastic scatter- 
ing of a photon off a charged object, and has proven to 
be very important as it provided early evidence that the 
electromagnetic wave is quantized [63, [68[. In QCD, the 
corresponding process is the gluon Compton scattering 
g + q(q) — > q(q) + 7*. The cross section in leading twist 
of pQCD [64| is given by: 



/ d 2 a(g + g) 
\dQ 2 d cosO 



GCS 



on-shell 



a 2 e 2 a s 



Q 2 



18Q 2 s 2 (l + cose) 
2s s - Q 2 



2Q 2 



2s 



(1 - cos 6) 



(1 + cos6) z 



(14) 



In terms of Mandelstam variables [20] it reads: 



/ d&(g + q) 
V dQ 2 dt 



GCS e 2 a 2 a S s 2 + t 2 + 2Q 2 u 



Jicll 



9Q 2 



-s 3 t 



xS(s + t + u- Q 2 ) , (15) 



which is obviously related by crossing transformation to 



III. DILEPTONS FROM DYNAMICAL 
QUASI-PARTICLES 

Let us now proceed to the calculation of the dilep- 
ton production by effective strongly interacting partonic 
quasiparticles with broad spectral functions. Dilepton ra- 
diation by the dynamical quasiparticles proceeds via the 
elementary processes illustrated in Figs.Q]and 2: the ba- 
sic Drcll-Yan q+q annihilation mechanism, Gluon Comp- 
ton scattering (q + g — > 7* + q and q + g — > 7* + q) , quark 
+ anti-quark annihilation with gluon bremsstrahlung 
in the final state (q + q — > .9 + 7*) and gluon decay 

(g-> q + q + l + 
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FIG. 3: (color online) Dimuon production cross sections in the Drell-Yan channel (q + q — > /i + +/i~). L.h.s. The cross section is 
presented versus the mass of the muon pair Q. The short dashes (black) line shows the on-shell, i.e. the standard perturbative 
result. The other lines show the off-shell cross section, in which the annihilating quark and antiquark have finite masses mi 
and mi with different values: mi = 0.3 GeV, mi — 0.15 GeV (solid magenta line), mi = 0.3 GeV, m,2 = 0.3 GeV (dash-dotted 
red line), mi = 0.3 GeV, mi = 0.6 GeV (dashed blue line), mi = 0.6 GeV, m-2 — 0.6 GeV (dash-dot-dot green line). R.h.s. 
The ratio of the off-shell cross section to the on-shell result for the different values of quark and antiquark masses is plotted 
versus Q — Qo, where Qo is the threshold value for the lepton pair mass. Line coding as in the figure on the l.h.s. 



A. Off-shell q + q in the Drell-Yan mechanism 

For the diagram (a) in Fig. [T] the off-shell cross section 
was obtained in [53j as: 



d 3 a(mi,m 2 ,pi,P2) 
dQ 2 dxpdq^ 



DY 



o f f — shell 



2Q 4 
xd(Q 2 

xS I xf 



Q (m 1 — Qm\m 2 + m 2 ) 



ml - ml- 2(pi -p 2 )) 



with 



2ira 2 e 2 



3Q 4 N c 8^( Pl -p 2 ) 2 -m\m 



(16) 



(17) 



In (|16p. iV c is the number of colors, e q (—e q ) is the frac- 
tional charge of the quark (antiquark), m 2 are the vir- 
tualities of the annihilating quark and antiquark, pi are 
their 4-momcnta. 

Just as in the standard pQCD 'on-shell' case, the mass 
of the produced dilepton pair is fixed to the invariant 
energy of the quark-antiquark collision: 

s = Q 2 . (18) 

However, the kinematical limit for the minimal dilepton 
mass is now higher than in the on-shell case: 



where mi ept is the mass of an electron or muon. Also, 
the incident current also changes in the off-shell case: 



J 



{ki-k 2 ) 2 — mfm 2 = — y (s — m 2 — m^Y—Am 2 ! 



s = (mi +m 2 + 2mi ept ) 2 > Am 2 ept , 



(19) 



compared to J = s/2 in the on-shell approximation. 

The approximation mi = m 2 — > in formula (|16|) is 
equivalent to restricting oneself to the leading term in 
the twist expansion, that is, in the case of the unpolar- 
ized Drell-Yan process, an expansion in powers of 1/Q. 
One can see that in this limit and additionally using the 
collinear kinematics pi± = p 2 j_ = we recover the stan- 
dard pQCD result Q. 

The off-shell cross sections are compared to the leading 
twist results in Fig. [3] Dimuon production cross sections 
in the Drell-Yan mechanism is plotted on the l.h.s. of 
Fig. [3] versus the mass of the muon pair Q = y/s. The 
short dashes (black) line shows the on-shell, i.e. the stan- 
dard perturbative result. The other lines show the off- 
shell cross section, in which the annihilating quark and 
antiquark have finite masses mi and m 2 with different 
values: m x = 0.3 GeV, m 2 = 0.15 GeV (solid magenta 
line), mi = 0.3 GeV, m 2 = 0.3 GeV (dash-dotted red 
line), mi = 0.3 GeV, m 2 = 0.6 GeV (dashed blue line), 
mi = 0.6 GeV, m 2 = 0.6 GeV (dash-dot-dot green line). 

The importance of higher twist corrections in the Drell- 
Yan process is illustrated by the ratio of the off-shell and 
on-shell integrated cross section a(Q) on the r.h.s. of 
Fig. [3] The ratio of the off-shell cross section to the on- 
shell result for the different values of quark and antiquark 
masses is plotted versus Q—Qo, where Qo is the threshold 
value for the lepton pair mass. Line coding as in the 
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figure on the previous figure. With increasing Q 2 , off- 
shell cross sections approach the leading twist - on-shell 
- result. 

In Ref. [60| , a model for implementing the higher twist 
corrections to the Drell-Yan process p + p — > l + + l~ + X 
was formulated making use of the above cross section, 
a Breit-Wigner parametrization for the quark spectral 
function and non-integrated parton distributions. The 
effect of quark off-shellness on the transverse momentum 
spectrum of the dileptons produced in pp collisions was 
found to be large at low and Q 2 . The calculations were 
compared to the data on the triple differential Drell-Yan 
cross section d 3 a/dQ 2 dxFdqr from experiment E866 [69| 
at Fermilab in pp collisions at 800 GeV incident energy. 
Both the slope and magnitude of the qr distribution of 
the Drell-Yan pairs were described by the adjustment of 
the spectral function parametrization. The distribution 
of the transverse momentum of lepton pairs produced in 
the Drell-Yan process off nuclei pA —> l + l~X was also 
reproduced within this model [53j j . However, the effect 
of quark off-shellness on the dilepton emission in heavy- 
ion collisions (A + A) has not been studied untill now. 
This question can be addresses by convoluting the off- 
shell cross section (|T6|) with effective spectral functions 
A(m,i) for quarks in plasma and with a model for the 
distribution of quarks in plasma (cf. Sections IV and V) . 



B. Off-shell gluon Bremsstrahlung qq — > gl + l 

In contrast to the off-shell cross section for the qq —> 
l + l~ process, the off-shell cross section for the qq anni- 
hilation with gluon Bremsstrahlung in the final state has 
not been calculated elsewhere. Therefore, we will provide 
here a short description of its evaluation. 

Starting from the formula for the unpolarized cross 
section 



da = 



EIAC/Peieall 



(2Tr) 4 5( Pl +p 2 -Z Pf ), (20) 



where the incoming quark and antiquark momenta are 
pi and p2, their masses mi and respectively; pf are 



the momenta of the outgoing particles, i.e. of the elec- 
tron (muon) and positron (anti-muon) and gluon, we note 
that the dilepton production cross section can be easily 
obtained from the cross section for the production of vir- 
tual photons as 



da(r) 



with 



do-(l+r) _ 
dQ 2 dt ~ 3ttQ 2 dt 



FF{Q 2 ,Ql) 



FF(Q 2 ,Q 2 ) (21) 




^5 ) . (22) 



where Q 2 = ^rn 2 ept , mi ept is the lepton mass. 

Furthermore, we define the momenta of the internal 
quark exchanged in the two relevant diagrams (see Fig. 1) 
as P3 = pi — q, p3 = pi — p2 — P3 and its mass as m 3 . The 
final gluon momentum is k and its mass is //. Then the 
matrix element of the process q + q — > g + 7* is 



M = M a + M b 



(23) 



where 
M a 



-e g eg s Tl 



e v {q)e a i(k) 



and 



Mb 



xui(pi,mi) [7"(p 3 + m 3 h a ]v :j (p2,m2) (24) 



T i e a i(k)e v (q) 
XUi(pi,mi) h v (p% + m 3 )Y] v j (p2,m 2 ),{25) 



e is the electron charge; e q is the quark fractional charge; 
Tjj is the generator of the SU(3) color group (that will 
yield the color factor in the cross section); e v {q) is the 
polarization vector for the virtual photon with momen- 
tum g; e a i{k) is the polarization vector for the gluon of 
momentum k and color /; Ui(p,m) is a Dirac spinor for 
the quark with momentum p, mass m and color i; and 
Vi(j>, m) is the spinor for the anti-quark. 



J 



The squared - and summed over spin polarizations as well as over color degrees of freedom 
be decomposed in the following summands: 



matrix element can 



]T \m\ 2 - Yl M a M * + E M b M » + E M * M >> + E M b 



M h 



(26) 



where the star denotes the complex conjugation. 

The spinors for quark states with mass m, contribute to the expression for the average matrix element only in the 
combinations ^2 u{p, TTii)u(p, mi) = (p+ TTii) (cf [70l |) and the correlation functions between the states with different 
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masses does not enter |M| 2 . Thus we find: 



2 q e 2 g 2 Tv{T 2 } 



01 - m l)(P3 2 - ml) 



[Tr {(p 2 - m 2 )la{pz + m 5 )^f u (pi + m 1 )-f' I {p 3 + m 3 )j"} 



1 

Q~ 2 

A 



Tr {(p 2 - m 2 )jcr(p~3 + m 3 )q(pi + mi)j a (p 3 + m 3 )q) 



k 2 



Tr 



.4 



k 2 Q 



|(p 2 - m 2 )k(p 3 + 1713)7* (pi, + mi)k(p 3 + m 3 )j^ 
Tr{(p 2 — m,2)k{p% + m 3 )q(p\ + mi)k(p 3 + mz)cj\ 



(27) 



_ . e 2 e 2 .g 2 Tr{T 2 } 

\ M a\ = 7Z3 Zl2TT~ I Tr {7<t(P3 + m 3 )j u (p\ + (p 3 + m 3 )^ a ' (p 2 ~ m 2 )} 



(pI ~ m 2 3 ) 2 



1 

Q 

A 



T^Tr {ja(p3 + m 3 )q(p\ + mi)q(p 3 + 7713)7° (p 2 - m 2 )} 



^Tr{fc(p 3 + m 3 )-f v {p\ + mi)Y{pa + m 3 )k(p 2 - m 2 )j 

Ac* - i 

^2^2 Tr |fc(p 3 + m 3 )q(p\ + m,i)q(p 3 + m 3 )k(p 2 - m 2 ) j 



(28) 



Note that by the transformation {p 3 —> p~ 3 ,pi — > p 2 ,p 2 —t Pi, Tlx — > — m 2 , m 2 —} —mi} we readily obtain ^ M£M a 
from M*M b and ^ |M fe | 2 from J] \MJ 2 ■ In equations (gTJ) and (J55J), A sets the gauge. We used the feynpar.m [71 1 
package of the Mathematica program [72j to evaluate the traces of the products of the gamma matrices. The resulting 
cross section is (here shown at A = 1): 



da(q q — > g l + l ) 



gBr 



a 2 a s e 2 8 (- Q 2 + s + t - /j 2 ) 5(s + t + u - m\ - m\ - Q 2 - u 2 ) 

dQ 2 dt J offsheU 27 Qi ^J_ 2 (mi 2 m2 2) + ( mi 2 + m2 2 _ s )2 s _ t j _ ^ ^ 2 

m 3 6 (si (s + t) + (s + 2t) 2 u+(s + 4t) u 2 -Q 4 (t + tt) + Q 2 (t 2 +u 2 )) 
+t u (-2 Q 6 (t + u) + 2 t u (2 s (s + t) + (2 s + t) u) + 4 Q 4 (t 2 + u 2 + s (t + u)) 
-Q 2 (2s 2 (t + u) + 4s (t 2 + u 2 ) + (t + u) (2t 2 -3tu + 2u 2 ))) +m 3 2 (2Q 6 (t + u) 2 

-Q 4 (t + u) (it (s + t) + (4s + t) u + 4u 2 ) + 2Q 2 {t 4 + t 3 u - t 2 u 2 + tu 3 + u 4 + s 2 (t + uf 

+2s (t + u) (t 2 +u 2 )) +tu (-3s 2 (t + u) - 2s (t 2 + 3iu + u 2 ) + 2 (t 3 + u 3 ))) - m 3 4 (- (s 2 (i - m) 2 ) 

+2 Q 6 (< + u) - 2 s (t + u) (t 2 -Ztu + u 2 ) +tu (3t 2 + 8tu + 3u 2 ) - Q 4 (5 i 2 + 2 1 u + 5 u 2 + 4 s (t + u)) 
+Q 2 (2 s 2 (t + u) + 4 s (t 2 + u 2 ) + 3 (t 3 + m 3 ))) - mi m 3 (< u (2 Q 4 {t + u) + s 2 (t + u) 

-Q 2 (t 2 + I2tu + u 2 + 3s (< + u))) +m 3 2 (t + u) (-2Q 4 (t + u) - s 2 (t + u) + Q 2 (t 2 + 12tu + u 2 

+3 s (t + u))) + m 3 4 (2 Q 4 (t + u) + s 2 (t + u) - Q 2 (3 s (t + u) + 7 (t 2 + u 2 )))) 
+ (-(m 3 6 (6Q 4 + 2s 2 + t 2 + u 2 -5Q 2 (t + u) + 4s (t + u))) 

+t u (-8 Q 6 + t u (2 s + t + u) + Q 4 (-4 s + 7 (t + u)) + Q 2 (4 s 2 + 2 (t - u) 2 + 5 s (t + u)) ) 

+m 3 4 (-8Q 6 + 2t(s + t) 2 + (2s 2 + 10st + t 2 ) u+(4s + t) u 2 + 2u 3 + Q 4 (-4s + 13 (t + u)) 

+Q 2 (4s 2 + 5s (t + u)-2 (t 2 + 8tu + u 2 ))) +m 3 2 (8Q 6 (t + u) + Q 4 (-It 2 - 20t u - 7u 2 + 4s (t + u)) 
-tu (2 s 2 + 6s (t + u) + 3 (t 2 + u 2 )) - Q 2 (4 s 2 (t + u) + (t + u) (4 1 2 - 13 1 u + 4 u 2 ) + s (7t 2 + 6tu + 7u 2 ))) 
+mim 3 (- (m 3 4 (lOQ 4 - 2 s 2 + t 2 + u 2 - s (t + u) + 2 Q 2 (2 s + t + u))) - t u (10 Q 4 - 2 s 2 + t 2 + u 2 
-s (t + u) + 2Q 2 (2s + t + u)) + m 3 2 (10 Q 4 (t + u) + (t + u) (-2 s 2 + t 2 + u 2 - s (t + u)) 
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+4Q 2 (-t 2 + 4tu-u 2 + s (t + u))))) u 2 + (m 3 2 -t) (m 3 2 - u) (4 mi m 3 Q 2 - 12 m 3 2 Q 2 - 10 Q 4 
-4 mi m 3 s + 4 m 3 2 s - 6 Q 2 s + 4 m 2 t + 3 Q 2 t - s t - t 2 + 4 to 3 2 u + 3 Q 2 u - s u - 4 1 u - u 2 ) u 4 
+ (m 3 2 - t) (m 3 2 - u) (2 (mi - m 3 ) m 3 - 2 Q 2 + t + u) /i 6 + 2 m 2 6 (m 3 2 - t) (m 3 2 - «) (2 Q 2 - s + u 2 ) 
+m 2 4 (to 3 2 - t) (m 3 2 - u) (2 mi m 3 Q 2 + 2 m 2 Q 2 + 4Q 4 - 6 Q 2 s + 2 s 2 - 5 Q 2 t + 2 at - 5 Q 2 u + 2 su 
-2 (m 3 (mi + m 3 ) + t + u) u 2 - 4/i 4 ) - 2m 2 5 (m 3 2 - t) (m 3 2 - u) (mi s + m 3 (3Q 2 - 2s + 3^ 2 )) 
-2 m 2 3 (m 3 2 - t) (m 3 2 - u) (- (mi s (-Q 2 + s + t + u - u 2 )) + m 3 (Q 4 

-3Q 2 {s + t + u) + 2s (s + t + u) -3 (s + t +it) /i 2 + ^i 4 )) - m 2 2 (m 3 2 - t) (m 3 2 - u) (2stu+ (4s 2 
-2 1 u + 3 s (* + u)) Ai 2 - 5 (2 s + t + u) u 4 + 6 ^ 6 + 5 Q 4 (i + u - 2 u 2 ) - 2 m 3 2 (Q - fi) (Q + /i) (2 Q 2 
-2s -t- u + 2/x 2 ) - Q 2 (5t 2 -2tu + 5u 2 + 5s (i + u) + 2 s /i 2 - 8 (t + u) // 2 + 2^ 4 ) +2mim 3 (s 2 

-(3s + i + w) /i 2 + 2/i 4 + Q 2 (-s + i + u-4^ 2 ))) -ma L (-2m 3 6 [q 4 -2Q 2 (s-5^ 2 ) + (s-^ 2 ) 2 ) 

+f it (2 s t u + ((i - it) 2 + s (f + it) J /i 2 - (t + u) fi 4 - Q 4 (t + u- A /i 2 ) + Q 2 ((t - uf + s (t + u) 

-26 (t + u) u 2 + 4u 4 )) + m 3 4 (2 stu + 2 s 2 (t + u) + (t - u) 2 u 2 - 3 s (t + u) fi 2 + (t + u) u 4 

+Q 4 (t + u + 4 u 2 ) + Q 2 (ft - u) 2 - 3 s (t + u) - 6 (t + u) u 2 + 4 u 4 ^j \ + m 3 2 (-2stu(s + t + u) 

-(t - uf (s + t + u) fi 2 + (t 2 + u 2 ) p 4 + Q 4 (t 2 + u 2 -4(t + u) fi 2 ) + Q 2 (- ((t - uj 2 (s + t + u)) 

-2 (5t 2 -52tu + 5u 2 ) /i 2 -4 (t + u) /it 4 ))) + m 3 (tu (-4Q 6 - 4stu + Q 4 (8 s + t + u) + Q 2 (~4s 2 + t 2 

+4tu + u 2 -s (t + u))) - (l2t 2 (s + t) 2 +t (-4Q 4 + 4s 2 + 19si+17i 2 + Q 2 (-I2s + 8t)) u+(l2s 2 

+8 Q 2 t + 19 s t - 4 1 2 ) u 2 + (24 s + 17 i) it 3 + 12 u 4 ) fi 2 + tu (4Q 2 + 8 s + t + u) p 4 -AtuyP 
+m 3 2 (-((t + u) (-4Q 6 -4stu + Q 4 (8s + t + u) + Q 2 (-4 s 2 + t 2 + Atu + u 2 - s (t + u)))) 
+ (-4Q 4 (< + u) + 28s 2 (t + u) + (t + u) (23t 2 + 8tu + 23u 2 ) +s (49t 2 + 74tu + 49u 2 ) 
+4Q 2 (t + u) (-3s + 2 (t + u))) u 2 - (t + u) (4Q 2 + 8s + t + u) u 4 + 4 (t + u) u 6 ) + m 3 4 (-4Q 6 
-4stu- (28s 2 + lH 2 + 32iit+llw 2 + 43s (t + u)) u 2 + (8s + t + u) u 4 - 4 u 6 + Q 4 (8s + t + u + 4u 2 ) 
+Q 2 (-4s 2 + t 2 +4tu + u 2 - 8 (t + u) u 2 + 4u 4 - s (t + u - 12 u 2 )))))] (29) 



In the off-shell case, the dilepton mass range is 4mf ept < 
Q 2 < (y/s — m q ) 2 , while the kinematical limits for the 
momentum transfer t are 

C?„ X = -f(Ci±C 2 ), (30) 

where 

Ci = 1 - (A + p 2 + ft + ft) + (ft - ft)(ft - ft), 

C 2 = y/(l - ft - ft) 2 - 4ftft 

x V(l - ft - ft) 2 - 4ftft (31) 

with 

ft = m 2 /s, ft = m 2 /s, ft = Q 2 /s, ft = u 2 /s. (32) 

Additionally, we note that there is a threshold in the 
CMS energy ^/s for the q + q interaction: 

s > max{(mi + m 2 ) 2 , (u + Q) 2 } . (33) 

One can easily check that the expression (|2"9"|) for 
rrii — > approaches the leading twist pQCD result, where 



I 

Mcwt = Z- 4 - We illustrate this in Fig. @] In Fig. 21 
the off-shell cross sections for the quark annihilation 
with gluon bremsstrahlung mechanism at various val- 
ues of quark and gluon off-shellnesses (masses) are com- 
pared to the on-shell (pQCD) result. Dashed black line 
shows the on-shell cross section with fj, cut = 0.206 GeV, 
red solid line presents the off-shell cross section for the 
gluon mass fixed to /i — 0.8 GeV and on-shell quark 
and anti-quark (mi = m 2 — m 3 = 0). Blue dash- 
dotted line gives the off-shell result for fi = 0.8 GeV, 
mi = m 2 = m 3 — m q = 0.6 GeV. One readily notices 
the shift of the maximum allowed mass of the pair to a 
lower value (in order to produce a massive gluon in the 
final state). For the rest of the Q values, the effect of 
the quark and gluon off-shcllncss reaches at most 50%, 
as is seen in the ratios of the cross sections, plotted in 
the r.h.s. of Fig.|H 

Next we compare the double differential off-shell and 
on-shell cross sections. The results for the transverse 
momentum distributions of the dileptons are presented 
in Fig. [5] Solid black line shows the differential on-shell 
cross section with \x cut = 0.206 GeV, blue dashed line 



10 




Q [GeV] Q [GeV] 

(a)Cross sections (b)Ratios of the cross sections 



FIG. 4: (color online) Comparison of off-shell and on-shell gluon Bremsstrahlung q+q — s> g+fi + )i~ cross sections at y/s = 4 GeV. 
L.h.s. Dashed black line shows the on-shell cross section with fj, cu t — 0.206 GeV, blue dash-dotted line presents the off-shell 
cross section for the gluon mass fixed to fj, = 0.8 GeV and on-shell quark and anti-quark (mi = mi = mz — 0). Red solid line 
gives the off-shell result for [s, = 0.8 GeV, mi = mi = mz = m q = 0.6 GeV. R.h.s. The ratio of off-shell to on-shell cross 
sections. Line coding is as on the l.h.s. plot. 
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Gluon Bremsstrahlung q+qbar->g+n, \i 
on-shell vs off-shell 




On-shell 

Off-shell, n=0.8 GeV, m =0 

Off-shell, ji=0.8 GeV, m =0.6 GeV 
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Gluon Bremsstrahlung q4-qbar->g+u |x 
on-shell vs off-shell 



• n=0.8 GeV, m =0.6 GeV 
q 

H=0.8 GeV, ra =0 



10 15 
q T 2 [GeV 2 ] 



20 



25 



30 



(a)Cross sections 



(b)Ratios of the cross sections 



FIG. 5: (color online) Comparison of the transverse momentum distributions of muon pairs produced in the gluon 
Bremsstrahlung q + q — > g + [i + fi~ channel in the off-shell and on-shell cases. L.h.s. Dashed black line shows the differ- 
ential on-shell cross section with fi cu t = 0.206 GeV, blue dash-dotted line presents the off-shell cross section for the gluon mass 
fixed to /i = 0.8 GeV and on-shell quark and anti-quark (mi = m,2 = mz = 0). Red solid line gives the off-shell result for 
fi = 0.8 GeV, mi = mi = 7713 = m q — 0.6 GeV. R.h.s. The ratio of off-shell to on-shell cross sections. Line coding is as on the 
l.h.s. plot. 



presents the off-shell cross section for the gluon mass 
fixed to [i = 0.8 GeV and on-shell quark and anti-quark 
(mi = ?Ti2 = rrij, = 0). Red dash-dotted line gives the 
off-shell result for /1 — 0.8 GeV, mi = 7712 = = m q — 
0.6 GeV. Again, we find the largest effect on the edge of 
the phase space, at the minimal qT- 



C. Off-shell gluon Compton scattering gq — > ql + l 

The cross section for the gluon Compton scattering can 
be calculated analogously to the calculation of the gluon 
Bremststrahlung cross section ((!?§]) above. On the other 
hand, the cross sections for g + q —> q + l + + l~ and 
g + q—i'q + l + + l~ is readily obtained from by the 
crossing transformation. 

Kinematic limits on s, t, Q 2 in the off-shell GCS pro- 
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cess are analogous to the q + q case. In the off-shell case, 
the dilepton mass range is 4m^ p( < Q 2 < (^/s — m q ) 2 , 
while the kinematical limits on the momentum transfer t 
are given by the formulae (|30[) - ([5Tj) with 

ft = mf/s, ft = fi 2 /s, ft = Q 2 /s, ft = m 2 2 /s, (34) 



while 



s > max{(mi + /i) 2 , (m 2 + Q) 2 }- 



D. Virtual gluon decay g — > qql + l 
quark decay q — > gql + l~ 



and virtual 



Although the process of real gluon decay g — > q + q + 
l = l~ is forbidden kinematically for perturbative parti- 
cles, it has a finite region of validity, if the gluon is off- 
shell, due to its broad with and finite pole mass. Anal- 
ogously, the virtual quark decay is also possible in the 
off-shell case. We present the Feynman diagrams for the 
corresponding processes in Fig. [51 The off-shell cross sec- 
tion for these processes can straightforwardly be obtained 
from (|29p by the crossing relation. For example, the cross 
section for q — > gql + l~ is obtained from q^gl + l~ by 
changing p 2 -» ~p 2 ■ 

Note that, in order to obtain the dilepton rates, the 
elementary cross sections have to be consequently con- 
voluted with the effective spectral functions for quarks 
and gluons. The magnitude and shape of the contribu- 
tions of the virtual decays to the dilepton rates is very 
sensitive to the final choice of the spectral function. In 
the DQPM the contribution from gluon decay is higher 
than that from virtual quark decay, since the gluonic 
quasi-particles are more massive and broader than the 
quarks [43j], and thus the kinematically allowed region 
is larger for the virtual gluon than for the virtual quark 
decay. 

On the other hand, within the DQPM parametriza- 
tions for the partonic spectral functions, both processes 
presented by the diagrams in Fig. [5] generate little dilep- 
ton yield anywhere except the extremely low masses: 
Q « 2m muon . Therefore, we refrain from plotting here 
the contributions explicitly and also do not consider them 
in the next section dedicated to the comparison of the 
yields from different mechanisms. 



IV. CONTRIBUTIONS OF THE DIFFERENT 
PROCESSES TO THE DILEPTON RATES FROM 

QGP 
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Model multiplicities 
for partonic collisions 



q+qbar 

g+q(qbar) 

Landau ansatz 
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FIG. 6: (color online) Multiplicities of q + q and g + q(g + q) 
collisions. The red solid line represents the multiplicities of 
the q + q collisions and the blue dashed line represents the 
probability for the gluon+quark(antiquark) interaction. The 
black dashed line shows for comparison the prediction in the 
Landau model for the heavy-ion collision - the power law fall 
as (v 7 *) -7 - 



n 1 r 



Dilepton rates 
on-shell, « s =0.4 



q+qbar->|x yi 

q+qbar->g+nn 

— q+g->q+nV 




FIG. 7: (color online) Dimuon rates dN/dQ from QGP calcu- 
lated using the cross sections in the on-shell approximation, 
as ~ 0.4. Black solid line shows the contribution of the Drell- 
Yan channel (q + q —¥ (i + fj,~), red dash-dotted line represents 
the contribution of the channel q+q — > g + blue dashed 

line shows the contribution of the channel q + g — > q + . 
The rates are in 'arbitrary units', reflecting our use of sim- 
plistic quark and gluon distributions in the QGP. 



In the following wc arc interested in the relative impor- 
tance of the different processes and their contributions 
to the dilepton yield of the strongly coupled quark-gluon 
plasma. 

Due to the factorization property proven in (22j . the 
dilepton emission from the QGP created in the heavy-ion 



collision is given by the convolution of the elementary 
sub-process cross sections (describing quark/gluon inter- 
actions resulting in the emission of dileptons) with the 
structure functions that characterize the properties and 
evolution of the plasma (encoded in the distribution of 
the quarks and gluons with different momenta and vir- 
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inverting the QCD EoS 
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1QCD 
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T/T 



FIG. 8: (color online) QCD energy density versus temper- 
ature from lattice QCD (sq uare symbols) [73|, the DQPM 
approach (red solid line) [43J and the 'AB-Bag' model (blue 
dashed line) [74| . The grey dash-dotted line shows the Stefan- 
Boltzmann limit. The arrow indicates the critical tempera- 
ture. 



tualities) 



d 3 CT QGP 



dQ 2 dxpdq 2 -, 

//>oo />oo />oo 

ds / dm\ / dm 2 / d\J F ab (s, m\, m 2 ) 
Jo Jo Jo 



xA ^ f ^ d 3 a abc (s,m\,m 2 ,n f ) 
dQ 2 dxpdq^ 



(35) 



where m\ and m 2 are the off-shcllnesses (i.e. virtuali- 
ties) of the incoming partons, ^ in the off-shellness of 
the outgoing parton, the indices a, b, c denote quark, an- 
tiquark or gluon so that all the considered mechanisms 
are covered. The cross sections a abc (s,m l1 , m l2 , /j,*) for 
the different processes have been derived in the previous 
Section. 

In (|35| . we integrate over the motion of partons, but 
also over their virtualities by employing phenomenologi- 
cal structure functions F ab that depend on the invariant 
energy s of the partonic sub-process as well as on the vir- 
tualities of the incoming partons and the spectral func- 
tion A(/^) for the virtuality distribution of the parton in 
the final state. Here should in principle be a two parti- 
cle correlator, but we work in the 2PI approximation so 
that the parton in the sQGP can be characterized by a 
single-particle distribution, therefore we assume that the 
plasma structure function can be approximated by 



F ab (s,mi,m 2 ) = A a (mi)A b (m 2 ) 



dNab 
ds 



(36) 



In this context, the quantity dN qq /ds has the meaning of 
the average multiplicity of q + q collisions in the plasma 
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FIG. 9: (color online) Temperature as a function of the energy 
density from lattice QCD (square symbols) [72], the DQPM 
approach (red solid line) [43] an d the 'AB-Bag' model (blue 
dashed line) (75 ]. 



as a function of the invariant energy of these collisions. 
Analogously, dN gq /ds denotes the multiplicity of g + q 
collisions. 

In order to address the relative importance of the dif- 
ferent mechanisms for the dilepton production in the 
sQGP, we need a quantitative model for the multiplici- 
ties of the collisions q + q and g + q (<? + <?) as functions of 
the center-of-mass energy yfs of these collisions. For this 
aim, we use the parametrizations, presented in Fig. [6] 
Inspired by the Landau model for heavy ion collisions, 
the parametrizations fall off with temperature according 
to the power low 0(T) ~ T~ r [H, while the magnitude 
of the distributions depend on the total available energy 
(~ snn)- Additionally, we notice that the energy of par- 
ton collisions is on average proportional to the tempera- 
ture y/s ~ T. In this admittedly very simplified picture, 
the functional form for dN/d^/s is taken as 



dN 
dyfs 



-iV oS 1/4 [(yCs-Pf + W] 



-3.5 



(37) 



By choosing (P qq = 0.5, W qq — 1.2) for qq collisions and 
(P gq = 1.2, W gq = 0.6) for gq collisions, we adjust the 
maximum of dN(q + q) to yfs ~ 0.5 GeV and shift the 
maximum of dN(g + q) to yfs « 1.2 GeV, reflecting the 
fact that the threshold y/so = m a + m b is higher for gq 
than for qq collisions. Indeed, the gluonic quasi-particles 
are expected to be more massive than the quark ones [43| . 
The parameter iVo remains to be fixed to experimental 
data or microscopic calculations, therefore we currently 
consider the distributions in "arbitrary units" . At high 
yfs the distributions approach the Landau ansatz const • 
T~ 7 , which is shown in Fig. [6] by the black dotted line. 

In Fig. [71 we plot the dilepton spectrum, assuming 
as = 0.4[77(. One observes that after the convolution 
with the distribution of possible yfs for the q + q annihi- 
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FIG. 10: (color online) Non-perturbative running coupling as 
a function of the local energy density e. The shadowed area in- 
dicates the energy densities reached in heavy ion collisions at 
SPS and RHIC. The arrow shows the critical energy-density. 



lation in sQGP, the yield of lepton pairs produced in the 
Bremsstrahlung process is below the leading order DY 
mechanism contribution. Thus, the DY rate is higher 
in the magnitude than that of the gluon Bremsstrahlung 
process despite the fact that the former one contributes 
only to the lepton pairs with a mass equal to the y/s. 

On the other hand, one notices from the Fig.[7]that the 
GCS mechanism is sub-leading, unless the gluonic con- 
tent of the plasma is orders of magnitude above the quark 
content, which is achieved neither at SPS nor at RHIC 
energies. A very high gluon content might be found at 
LHC, in which case, the GCS process would give a con- 
siderable contribution to the dilepton yield of the QGP. 

Let us remind that the running coupling as depends 
on the local energy density e. The DQPM [52| provides a 
good parametrization of the QCD running coupling as a 
function of temperature in the non-perturbative regime 
for temperatures close to T c (cf. Fig. 1 in Ref. [52|), Note 
that close to T c the full coupling calculated on the lattice 
increases with the decreasing temperature much faster 
than the pQCD prediction. 

The relation between the energy density and temper- 
ature, i.e. the equation of state, has also been extracted 
on the lattice [fg] - see Fig. H The DQMP model de- 
scribes the QCD equation of state even at T ~ T c [43j . 
A rather simple parametrization for the QCD equation 
of state - the "AB-Bag model" - is proposed in 74| and 
provides a good fit of the SU(3) lattice data above T c ; 
we extend this model to 3-flavors and also compare to 
lattice data [73| in Fig [5] 

In Fig [U] we reversed the relation and present the tem- 
perature as a function of the energy density. Using this 
relation, we obtain the running coupling as a function 
of the energy density e; we present as vs. e in Fig. [TO] 
On the other hand, simulations in transport theory |75j 
have shown that the local energy densities achieved in the 
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FIG. 11: (color online) Dimuon rates dN/dQ from QGP cal- 
culated using the cross sections in the on-shell approximation, 
as = 0.8. Black solid line shows the contribution of the Drell- 
Yan channel (q + q — > fi + fj,~), red dash-dotted line represents 
the contribution of the channel q + q — ¥ g+fJ,~ M + , blue dashed 
line shows the contribution of the channel q + g — > q + . 
The rates are in 'arbitrary units', reflecting our use of sim- 
plistic quark and gluon distributions in the QGP. 



course of heavy-ion collisions at SPS and RHIC energies 
reach at most 20 GeV/fm 3 ; this region is high-lighted in 
Fig. [10] by a shadowed area. 

One observes that as at the energy densities of interest 
is on average approximately 0.8. Using this value for as, 
we compare the rates in Fig. [TTJ In this case, the contri- 
bution of the 0(as) diagrams (gluon-Compton scatter- 
ing qg — > qj* and gluon Bremsstrahlung qq — > gj*) is no 
more subleading to the DY annihilation mechanism! 

Next, we plot the dilepton rates - within our approxi- 
mate model for the parton collision density in the plasma 



for the case of massive quarks and gluons in Fig. 12(a) 



The rates are calculated by convoluting the off-shell cross 
sections obtained in the previous Section with model 
dN/d^/s. Quarks and gluons are massive quasi-particles, 
quark masses being set to m q = 0.3 GeV and gluon mass 
fi = 0.8 GeV. Black solid line shows the contribution 
of the Drell-Yan channel (q + q —> /i + /i~), red dash- 
dotted line represents the contribution of the channel 
q + q — > g + blue dashed line shows the contri- 

bution of the channel q + g — > q + fi + \T~ . The rates in 
the three channels are modified compared to the mass- 
less case (cf. Fig. [TTj) . In particular, we point to the 
clear threshold behavior of the Born (Drell-Yan) term, 
noted also by the authors that have studies the dilepton 
production by massive partons previously (43l - [45| . 
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FIG. 12: (color online) Dimuon rates from QGP dN/dy beyond the on-shell approximation, as = 0.8. L.h.s. dN/dQ calculated 
using the derived off-shell cross sections for quarks and gluons as massive quasi-particles, quark mass being set to m q — 0.3 GeV 
and gluon mass to /i = 0.8 GeV. Black solid line shows the contribution of the Drell-Yan channel (q + q — > /i + ^t~), red dash- 
dotted line represents the contribution of the channel q + q — > g + n~ n + , blue dashed line shows the contribution of the channel 
q + g — > q + fi + fi~ . R.h.s. dN/dQ calculated in the fully off-shell case of massive and broad dynamical quasi-particles. The 
rates are calculated using the derived off-shell cross sections and convoluted with effective spectral functions of the Breit-Wigner 
type. The parameters of the spectral functions are: the peak of the quark spectral function is located at 0.3 GeV, the width is 
r = 0.3 GeV, the peak of the gluon spectral function is at 0.8 GeV, the width to T = 0.3 GeV. Line coding is as on the l.h.s. 
The rates on both plots are in 'arbitrary units', reflecting our use of simplistic quark and gluon distributions in the QGP. 



V. EFFECT OF THE FINAL QUARK AND 
GLUON WIDTH ON THE QGP RADIATION 

Finally, we calculate the QGP dilepton rate, taking 
into account not only the finite masses of the partons, 
but also their broad spectral functions, i.e. finite widths. 
For this purpose, we convolute the off-shell cross sections 
obtained in section IIIII with dNj dyfs and with the spec- 
tral functions A(m) according to the equation (1351) . 

The partonic spectral functions are related to the imag- 
inary part of the trace of the effective propagator and 
to the partonic self-energies £ as follows: 



the width is T = 0.3 GeV, the peak of the gluon spectral 
function is at 0.8 GeV, the width T = 0.3 GeV.. The rate 
from the Drell-Yan mechanism is shown by solid black 
line, while for the Q(cts) 2 — » 2 processes by the dashed 
blue and dash-dotted red lines. We have checked nu- 
merically that the rates for different values of the gauge 
parameter lie on top of each other. For comparison, the 
rates shown on the l.h.s. of Fig. [T2]have been obtained in 
the previous Section in the approximation of zero width 
(6- functions for spectral functions). 



A( P ) = -QDftp) 

7T 



9E(p) 



7T [p- 



^current 



ReY,(p)] 2 + [9fE(p)] 



r(38) 



For the current qualitative analysis, we use the approxi- 
mation of constant real and imaginary parts of the self- 
energy, which corresponds to constant finite average mass 
for quarks (< m q >) and gluons (fi) and their width T. 
Within these approximations, the spectral function has 
the Breit-Wigner form. 

The results of the numeric convolution of the off- 
shell cross sections with the spectral functions and the 
dN/dy/s are shown in Fig. 12(b) for realistic values of 



the spectral function parameters inspired by DQPM: the 
peak of the quark spectral function is located at 0.3 GeV, 



By dressing the quark and gluon lines with effec- 
tive spectral functions we model the effect of the quasi- 
particle interaction, including their multiple scattering. 
One observes by direct comparison of the l.h.s. and r.h.s. 
of the Fig.[T2]that the effect on the dilepton rates is quite 
dramatic. In particular, the threshold of the Drell-Yan 
contribution is "washed-out" . In this observation we con- 
firm the results of [5l|. On the other hand, the effect of 
the partonic width and/or of multiple scattering on the 
2 — > 2 processes has not been studied so far. We ather th e 
predicted shape of the dilepton spectrum in Fig. 12(b) is 



realized remains to be answered in the comparison to ex- 
perimental data [76|. 
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VI. SUMMARY AND OUTLOOK 

In the present work, we have derived the off-shell cross 
sections of dilepton production by the constituents of the 
strongly interacting quark-gluon-plasma in the reactions 
5 + (Drell-Yan mechanism), q + q g + 

(quark annihilation with the gluon Bremsstrahlung in the 
final state), q(q) +g — > q(q) + l + l~ (gluon Compton scat- 
tering), g ->■ q + q + and q(q) q(q) + g + 
(virtual gluon decay, virtual quark decay) by dressing 
the quark and gluon lines in the perturbative diagrams 
with effective non-perturbative propagators. 

We show that finite quark and gluon virtualities mod- 
ify the magnitude, Q and qT dependence of the cross 
sections of all the calculated processes compared to the 
leading twist perturbative results for massless partons. 
The modification is higher at low Q 2 and on the edges of 
the phase space. 

In QGP, quarks and gluons quasiparticles have broad 
spectral functions (due to large interaction rates) , where 
the width encodes the effects of rescattering and gluon 
radiation. We have found that the finite width of partons 
has a dramatic effect on the dilepton rates in QGP, espe- 
cially at low Q 2 . In particular, the threshold of the Born 



term is smeared, while the contribution of the 2—^2 
processes to the rate is increased. 

The cross sections obtained in this study will form the 
basis of a consistent calculation of the dilepton produc- 
tion in heavy ion collisions at SPS and RHIC energy by 
implementing the partonic processes into a transport ap- 
proach of the PHSD collaboration (76|. The comparison 
to the dilepton data of the NA60 and PHENIX Collab- 
orations double differentially in mass and pt will open 
the possibility to study the relative importance of dif- 
ferent processes in the dilepton production and guide us 
towards a better understanding of the properties of mat- 
ter created in heavy-ion collisions. 
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